Title: "An Apparatus and method for monitoring a vascular access 
of a patient subjected to an extracorporeal blood treatment." 

DESCRIPTION 

Background of the Invention. 

5 The invention relates to a system for monitoring vascular access 
in a patient undergoing extracorporeal blood treatment. 

Specifically, though not exclusively, the invention can be 
usefully applied in the field of extracorporeal treatment for 
kidney failure. 

10 Setting up an extracorporeal blood treatment, such as. for 
example hemodialysis therapy, requires blood circulation in an 
extracorporeal circuit connected to the cardiovascular circuit 
of the patient through a vascular access . 

The blood, taken from the patient and sent through an 
15 extracorporeal circuit, is subjected to a treatment, generally 
passing through a treatment unit (for example a dialyzer filter) 
and thereafter returned to the patient. The vascular access is 
where the blood is removed from the cardio-vascular system of 
the patient and returned to the system. 

20 One of the vascular accesses most commonly used in hemodialysis 
therapy is the Cimino-Brescia artero-venous fistula access. 
Other vascular access types are known, however. For reasons of 
simplicity the present description will make reference to the 
artero-venous fistula as an example of vascular access, without 

25 excluding other types of vascular access from the claimed field 
of protection. 

In an extracorporeal treatment the blood is usually taken from 
the vascular access by an arterial needle fluidly connected to 
the extracorporeal circuit. After having passed through the 
30 treatment unit, the blood is sent back to the vascular access 
through a venous needle. Generally blood circulation in the 
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extracorporeal circuit is performed by a positive displacement 
pump, generally peristaltic. 

One of the problems of extracorporeal blood treatment is 
monitoring the efficiency of the vascular access. 

5 A parameter indicative of this efficiency is the blood flow rate 
which the vascular access can supply. This flow rate is usually 
greater than the blood flow rate through the extracorporeal 
circuit. For example, in normal conditions the blood flow rate 
at the vascular access is about 800-M200 ml/minute, while the 
10 blood flow rate in the extracorporeal circuit varies between 100. 
and 70 0 ml /minute. 

The flow rate at the vascular access can diminish due to a 
vascular pathology, such as, for example, a stenosis, i.e. a 
narrowing of the blood passage section, or for example due to a 
15 drop in cardiac output. The presence and location of a stenosis 
at the vascular access should be determined as soon as possible 
in order to prevent the stenosis degenerating into a thrombosis 
(occlusion of the blood vessels) . 

A reduced- ef f iciency vascular access can lead to the undesirable 
20 phenomenon of recirculation of' blood during treatment. 
Recirculation consists in the presence during treatment of blood 
flow proceeding in an opposite to the desired direction, i.e. 
from the return zone of the treated blood (venous needle) to the 
supply zone of the blood to be treated (arterial needle) . Flow 
25 recirculation therefore consists in return (in the 
extracorporeal circuit) of blood which has already been 
subjected to treatment, with a consequent diminution in 
treatment efficiency. 

Various systems have been proposed for monitoring vascular 
30 access and, more generally, the cardiovascular system of a 
patient subjected to extracorporeal blood treatment. 

EP 1 044 6 95 A2 teaches a method for determining the blood flow 
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rate in a vascular access during hemodialysis treatment. The 
method varies the blood flow rate of the extracorporeal circuit 
and measures the arterial and venous pressures in the 
extracorporeal circuit during the above-cited flow rate 
5 variations. The operations are carried out in two different 
conditions: first with the vascular access open, in which a part 
of the blood flow passes through the vascular access between the 
withdrawal needle and the return needle, and then when the 
vascular access is closed, in which the vascular access flow 

10 between withdrawal needle and return needle is zero. According 
to the method of EP 1 044 695 A2 , vascular access blood flow 
rate, with the vascular access open, is judged to be equal to 
the blood pump flow rate at which the difference of arterial 
pressure (or venous pressure) in the two different situations is 

15 zero. 

This method has the drawback that it is necessary to intervene 
mechanically on the fistula to interrupt blood flow. 

WO 00/18451 teaches a method for determining the flow in a 
fistula of a patient using an extracorporeal blood flow circuit, 

20 such as for example a hemodialysis circuit, in which the blood 
flows from a withdrawal point in the fistula to a return point 
in the fistula. The method varies the blood flow rate in the 

. extracorporeal circuit and takes a reading of a signal which can 
be correlated with the fistula flow rate downstream of the 

25 withdrawal point, The blood flow rate upstream of the withdrawal 
point is evaluated at equal the blood flow rate obtaining in the 
extracorporeal circuit when the fistula blood flow rate 
downstream, read with the above -described signal, is zero. WO 
00/18451 includes an embodiment in which the signal which can be 

30 correlated to the fistula blood flow downstream of the 
withdrawal point is generated by an ultrasonic sensor which 
operates directly on the patient's vascular access. 

The use of a sensor to directly measure the blood flow rate. in. 
the fistula tract comprised between the withdrawal needle and 
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the return needle leads to a certain constructional 
complication, as well as some discomfort for the patient. 

EP 1 020 199 A2 teaches . a method for detecting the presence of a 
stenosis in a vascular access during extracorporeal blood 
5 treatment . The method includes the use of at least one pressure 
sensor predisposed in the extracorporeal circuit along the 
arterial line upstream of the blood pump. A stenosis can be 
calculated from the entity of the pressure pulse measured by the 
pressure sensor. 

10 A pressure sensor can be placed on the arterial line too 
downstream of the blood pump and upstream of a dialyzer, and a 
further pressure sensor can be placed on the venous line 
downstream of the dialyzer. The method also includes a reading 
of the pressure pulse frequency and use of that frequency as a 

15 signal entity correction factor. The pressure pulse frequency 
signal can be corrected by means of a function depending on the 
blood pump flow rate. 

The data deductible from the method described in EP 1 020 199 A2 
is however limited: in particular, the method provides only a 
20 general indication of the hemodynamic state of the fistula, 
signalling the presence of a stenosis, but it cannot gather more 
detailed data, such as for example the vascular access blood 
flow rate or the location of any stenoses found. 

US 5,830,365 teaches a method for determining some hemodynamic 
25 parameters, among which the blood flow rate in a fistula during 
an extracorporeal blood treatment. The method involves the 
alteration of at least one chemical -physical characteristic of 
the blood in the venous line of the extracorporeal circuit, and 
recording the change which occurs in the arterial line following 
30 this alteration. The alteration can be a change in the 
concentration of an indicator, or a change in the temperature or 
pressure. In a specific embodiment use is made of a hemodialysis 
machine provided with a dialyzer where a dialysis solution 
containing an indicator flows and the concentration change of 
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the indicator in the venous and arterial lines of the 
extracorporeal circuit connected to the dialyzer is registered. 
In the venous line the concentration of the indicator increases 
by effect of back- filtration through the dialyzer. In the 
5 arterial line the concentration of the indicator increases by 
effect of recirculation in the fistula. The change of 
concentration in the arterial and venous lines is read by 
ultrasonic sensors. Alteration (in this case the change . in 
concentration) is performed in two stages: first when the blood 

10 flows in the normal, direction through the extracorporeal 
circuit, then when the blood flows in the opposite direction. 
The method includes the use of a device for inverting the blood 
flow direction in the extracorporeal circuit. According to the 
method taught in US 5,830,365 the change in concentration 

15 measured in the first stage enables calculation of recirculation 
at normal flow rate, while the change in concentration measured 
in the second stage enables calculation of recirculation when 
the flow is inverted. The two calculated values thus enable a 
calculation of various hemodynamic parameters among which the 

20 blood flow rate in the fistula. 

However the alteration of the chemical -physical properties of 
the blood and the inversion of the flow during the course of 
extracorporeal treatment lead to various drawbacks: a 
constructional complication, a delay in carrying out the 
25 treatment, an invasive intervention on the blood, quite removed 
from the course of normal treatment. 

WO 02/04044 teaches a method for identifying problems in 
arterial flow during an extracorporeal blood treatment in which 
the blood is transferred, by means of a positive displacement 

30 pump, from the vascular access of a patient to a blood treatment 
device through an arterial line and then sent by the treatment 
device to the vascular access through a venous line of the 
extracorporeal circuit. The method consists in measuring the 
amplitude of the periodic variations in pressure in the venous 

35 line induced by the rotation of the blood pump, by comparing the 
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variations with a threshold value and generating a control 
signal if the threshold value is exceeded. WO 02/04044. further 
describes another method according to which, during a dialysis 
treatment, the amplitude of the periodical variations of the 
5 pressure of the dialysis fluid (and not the venous line) is 
measured. The result is compared with a threshold value, and if 
the threshold value is exceeded a control signal is generated. 

The methods described in WO 02/04 044 are not however able to 
provide data relating to the blood flow rate at the vascular 
10 access. 

US 6,221,040 discloses a system for monitoring a vascular access 
during a dialysis treatment in which the pressures in both the 
arterial and venous branches of the extracorporeal blood system 
are monitored by piressure sensors. A computer unit generates 

15 first and second characteristic values for the integrity of the 
vascular access from measured arterial and venous pressures. An 
analyser unit analyses the integrity of the vascular access by 
comparing the first and second characteristic values to first 
and second ranges of predetermined values. Calculating a sum of 

20 the venous and arterial pressure generates the . first 
characteristic value, and calculating a difference between the 
venous and the arterial pressure generates the second 
characteristic value. 

The object of US 6,221,040 is to provide a monitoring system 
25 that allows detection of the venous cannula slipping out of the 
vascular access as well as detection of a blood leak in the 
venous branch of the extracorporeal circuit. It is not directed 
to determination of fistula flow. 

US 5,866,015 and EP 0 773 035 disclose a method for determining 
30 hemodynamic parameters during an extracorporeal hemotherapy, 
including the steps of measuring the blood temperature in the 
arterial branch of the extracorporeal circuit, varying the blood 
flow in the extracorporeal circuit, storing the values of the 
extracorporeal blood flow and the measured values of the blood 
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temperature, and determining a value of the blood flow from the 
stored sequence of value pairs of blood temperature, and of 
extracorporeal blood flow, at which value, after it is exceeded, 
the amount of the change in the blood temperature within a 
specific blood flow interval is greater than a predetermined 
limiting value. The fistula flow is inferred from the determined 
blood flow value. 

The method is based on the fact that, the measuring curve 
existing in discrete measured values is able to be represented 
by two subf unctions, the first subf unction indicating the blood 
temperature as a function of the extracorporeal blood flow for 
blood flow values smaller than the fistula flow or equal to the 
fistula flow, and the second subf unction indicating the blood 
temperature as a function of the blood flow for blood flow 
values greater than or equal to the fistula flow. The 
intersection of the two subfunctions indicates the point where 
the extracorporeal blood flow equals fistula flow. Thus, from 
the "break point" of the characteristic function curve, i.e., 
from the discontinuity in the rise of the curve, the point is 
able to be defined where fistula recirculation begins, i.e., 
where blood flow equals fistula flow. 

In addition to measuring temperature, the concentration of a 
blood constituent (hematocrit) can also be measured, as can the 
density, speed of sonic propagation, optical density, and 
conductivity or viscosity. 

The blood characteristic to be measured must have a different 
value in the venous branch of the extracorporeal circuit than it 
does in the blood flowing to fistula. It is assumed that the 
blood characteristic, preferably the temperature, is kept 
constant in the venous branch of the extracorporeal circuit 
while the measured values are recorded. If this characteristic 
is not constant, a regulating device to keep the characteristic 
in the venous branch constant must be provided. In the case of a 
temperature measurement, for example, this can be realized as a 
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temperature controller . 

Another drawback of this method is that the delivery rate of 
. , blood pump,, which predetermines the extracorporeal blood, flow, 
is increased starting, from a lower value to an upper limiting 
5 value which must be greater than the fistula flow to be 
expected. Fistula flows can only be determined within the 
adjustable blood flow range. Therefore the fistula flow is not 
determinable if it is equal to or greater than the upper 
limiting value of the adjustable blood flow range. ■ 

10 The prior art comprises the scientific publication entitled: 
"On-line dynamic measurement of fistula pressure during 
hemodialysis for detection of access stenosis and bad needle 
placement", Abstract from the 24th EDTNA-ERCA Conference, 
Prague, 5-8 July 1997, page 23, authors Polaschegg, Techert and 

15 Wizemann. 

According to this publication it is possible to calculate . the 
pressure of a vascular access by measuring the pressure in an 
extracorporeal blood circuit connected to the vascular access, 
with the aim of detecting any stenoses in the access itself. 

20 In a scientific publication entitled "Dynamic pressure 
measurement for detection of blood access stenosis", published 
in the EDTNA-ERCA Journal, 1998, XXIV, 4, on pages 3 9-44, 
authors Polaschegg, Techert and Wizemann, more detail is given 
on monitoring problems in a patient's vascular access. The 

25 method is based on the determination of the venous and arterial 
pressures (upstream of the blood pump) in an extracorporeal 
blood circuit connected to the vascular access to be monitored. 
The method comprises a preliminary stage in which, through in 
vitro tests in which the extracorporeal circuit is not connected 

30 to a real vascular access, fluid resistances in the arterial and 
venous lines of the extracorporeal circuit are calculated. 
During a second stage the extracorporeal circuit is connected to 
the real vascular access of the patient in order to initiate an 
extracorporeal treatment . During the extracorporeal treatment 
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the venous and arterial pressures are calculated in the 
extracorporeal circuit . As the venous and arterial pressures in 
the extracorporeal circuit are known, as are the fluid 
resistances in the arterial and venous lines of the 
5 extracorporeal circuit, the pressures in the vascular access can 
be calculated. The dynamic measurement at different flow rates 
and the comparison with static measures enables stenoses at the 
vascular access to be identified. 

The scientific publication entitled "Pressure drops in cannulas 
10 for hemodialysis", author H.D. Polaschegg, published in The 
International Journal of Artificial Organs, Vol. 24, No. 9, 

2001, pp. 614-623, relates to a method for determining a fall in 
pressure in an arterial or venous line in hemodialysis, with 
which the vascular access pressures can be determined starting 

15 from the pressures measured in the extracorporeal circuit of the 
hemodialysis machine . 

The scientific publication entitled "Extracorporeal pressure 
monitoring and the detection of vascular access stenosis", 
authors Kleinekof ort , Kraemer, Rode and Wizemann, published in 
20 The International Journal of Artificial Organs, Vol. 25, No. 1, 

2002, pp. 45-50, presents a method for identifying the presence 
of stenoses in a vascular access, even where the stenosis is 
located between the withdrawal needle and the return needle. The 
method comprises measuring the static pressures in the arterial 

25 and venous lines of an extracorporeal circuit and in calculating 
the pressures at the vascular access at the withdrawal point and 
the return point. These pressures, which correspond to the 
pressures which would be measured by two pressure sensors 
directly connected to the withdrawal and return needles, are 

30 used in order to identify the presence of a stenosis. A 
knowledge of the pressures both at the point of withdrawal and 
at the point of return of the vascular access provides more 
accurate indications and enables a first approximate 
localization of the stenosis, especially enabling to detect if 

35 the stenosis is in venous tract or is located between the 
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needles . 



The method described here is not however able to determine the 
.blood flow rate in the vascular access. 

The publication entitled "Utility of intra-access pressure 
monitoring in detecting and correcting venous outlet stenoses 
prior to thrombosis", in Kidney International, Vol. 47 (1995) , ; 
pages 13 64-13 73, authors Besarab, Sullivan, Ross, Moritz, 
teaches a method for deriving the pressure internally of the 
vascular access (intra-access pressure) from the pressure 
measured in the hemodialysis machine, as a function of the type 
of needle used, the blood flow rate of the hemodialysis machine, 
and the hematocrit of the blood. Other methods for determining 
the pressure at the vascular access are cited or described in 
the following publications: 

"Detection of access strictures and outlet stenoses in 
vascular accesses. Which test is best?", in ASAIO Journal, 
1997, Vol. 43: pages M543-M547, authors Besarab, Lubkowski, 
Frinak, Ramanathan, Escobar ; 

"Simplified measurement of intra-access pressure", in 
Journal of the American Society of Nephrology, 1998, Vol. 
9, pages 284-289, authors Besarab, Frinak, Sherman, 
Goldman, Dumpier, Devita, Kapoian, Al-Saghir, Lubkowski; 

"Effect of systemic hemodynamics on flow within vascular 
accesses used for hemodialysis", in ASAIO Journal 2001, 
Vol. 47, pages 501-506, authors Besarab, Lubkowski, Vu, 
Aslam, Frinak; 

"Dynamic venous access pressure ratio test for hemodialysis 
access monitoring", in American Journal of Kidney Disease, 
Vol. 40, N° 4 (October), pages 760-768, 2002, authors 
Frinak, Zasuwa, Dunfee, Besarab, Yee . 

An abstract entitled "A novel model -based method for monitoring 
the hemodialysis vascular access", published in the Journal of 
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the American Society of Nephrology, 2001, Vol. 12, N. A1513, 
pages 294A-295A, authors Lodi , Monari , Fava, Paolini, Grandi , 
Galato, Cavalcanti, cites a mathematical model based on the 
hemodynamic description of the vascular access which enables the 
5 arterial and. venous pressures at the vascular access to be 
calculated and also the flow in vascular access starting from 
extracorporeal arterial and venous pressures. The model, which 
includes three parameters (resistance to flow of the 
anastomosis, resistance between arterial and venous access, the 

10 resistance which expresses the efficiency of venous circulation 
drainage) , was used to analyse the data gathered during a normal 
hemodialysis therapy operation. The abstract states that the 
extracorporeal venous and arterial pressures were measured after 
having set four different flow rates on the blood pump and that 

15 the above-cited parameters included in the mathematical model 
were calculated using the mathematical model. 

Summary of the Invention. 

The present invention provides a system for controlling vascular 
access adequacy during an extracorporeal blood treatment. 

20 An aim of the invention is to enable calculation of some 
hemodynamic parameters at the vascular access. Knowledge of 
these parameters enables both regulation of the blood pump flow 
rate operation in the extracorporeal circuit and intervention in 
case of detection of a pathological situation in the vascular 

25 access. 

A further aim of the invention is to enable evaluation of the 
blood flow circulating in the vascular access of a patient 
during an extracorporeal blood treatment. 

A further aim of the invention is to make available a system for 
30 evaluating vascular hydraulic resistance in various tracts of 
the patient's vascular system. In particular, an aim of the 
invention is to evaluate vascular resistance upstream of the 
blood withdrawal zone from the vascular access, downstream of 
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the blood return zone, and in the tract of vascular access 
comprised between the withdrawal zone and the return zone. 

, An advantage of the invention is that it provides indicative 
values of the efficiency of the vascular access simply, 
5 automatically, using devices (such as for example pressure 
transducers, blood pump, drainage pump) which are normally 
already present in machines for extracorporeal blood treatment. 
A further advantage is that the invention enables monitoring of 
the vascular access at any time during the extracorporeal blood 
10 treatment . 

A further advantage of the invention is that the monitoring, 
procedure does not cause extra stress to the patient. The 
procedure can be carried out by means of variations in the blood 
pump or the drainage pump flow rates, or both, within flow rate 
15 intervals which are normally compatible with the extracorporeal 
treatment the patient undergoes. The intervals can be those 
normally used during the course of therapy. 

These aims and others besides are all attained by the invention 
as it is characterised in one or more of the appended claims. 

20 In a special function of the invention, a mathematical model is 
used which contains at least two parameters in which a first 
parameter relates to the hemodynamics of the vascular access, 
and a second parameter relates to the blood flow rate in the 
extracorporeal circuit . 

25 The mathematical model comprises a third parameter relating to 
at least one blood characteristic: this characteristic can be 
any physical, chemical or physical -chemical property thereof 
which characterises the blood in a vessel and which can be 
related to the blood flow rate in that vessel. A peculiarity of 

30 the invention is that the mathematical model used describes the 
relationship between the selected blood property (physical, 
chemical or physical - chemical ) and the blood flow rate in the 
vessel. In particular the mathematical model describes the 
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relationship in the vascular access. For example, the 
mathematical model can describe the fluid-dynamic situation of 
the vascular access; the model can describe a relationship 
between the difference of pressure at two points of the vascular 
5 access and the flow rate crossing the points. Apart from the 
pressure it is also possible to select other properties 
(physical, chemical or physical -chemical) of the blood which are 
. influenced by the flow rate, such as, for example: the. 
difference in induced potential, speed of sound, optical 
10 characteristics, temperature, concentration of an indicator, and 
so on. 

According to the invention, the monitoring of the vascular 
access is performed by varying the flow rate of at least one 
fluid (for example blood or the product of ultrafiltration) , 
15 which runs either in the extracorporeal circuit or in at least 
one hydraulic line (for example an ultrafiltration line) 
connected to the extracorporeal circuit. 

The monitoring can be carried out by varying both the above - 
cited flows. 

20 The monitoring determines the values of at least one 
characteristic of the blood, in at least one zone of the blood 
circulation path, and at at least two different values of the 
flow rate of the fluid. 

As mentioned above, the cited characteristic of the blood can be 
25 a physical, chemical or chemical -physical one. In an embodiment 
of the invention, among the various characteristics of the blood 
that depend on blood flow, the selected characteristic to be 
used is the pressure. 

The monitoring procedure involves calculating one or more of the 
30 hemodynamic parameters of the vascular access contained in the 
mathematical model, using the values of the blood characteristic 
determined previously during the course of the procedure. 
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In an embodiment of the invention, a multiplicity of values of 
the blood characteristic is determined; then the said 
hemodynamic parameters are calculated, by means of the 
mathematical model,. using approximation algorithms (of known 
5 type) . The algorithms can be chosen, for example, from those 
which enable determination of the value of the hemodynamic 
parameter, by virtue of which the blood characteristic . values 
calculated using the mathematical model, at different flow rate 
values, are those which are closest to the blood characteristic 
10 values which were previously determined during the. course of the 
procedure, at the same flow rate values. 

In an embodiment of the invention, the mathematical model used 
is descriptive of the pressure variation at the vascular access: 
it comprises at least one hemodynamic parameter relative to at 
15 least one characteristic of the vascular access; at least one 
parameter relative to the blood characteristic; and at least one 
parameter relative to the blood flow rate in the extracorporeal 
circuit . 

The hemodynamic parameter can be relative to at least one of the 
20 following characteristics of the vascular access: the blood flow 
rate upstream of a withdrawal zone of the blood from the access, 
the blood flow rate between the withdrawal zone and a blood 
return zone at the access, the blood flow rate downstream of the 
blood return zone, the vascular hydraulic resistance upstream of 
25 the blood withdrawal zone from the access, the vascular 
hydraulic resistance between the blood withdrawal zone and the 
blood return zone, and the vascular hydraulic resistance 
downstream of the blood return zone. 

In a further embodiment of the invention, the monitoring 
30 procedure includes determining the values assumed by the blood 
characteristic in at least two zones of the blood circulation 
path (where the blood circulation path comprises both the 
intracorporeal circuit and the extracorporeal circuit) and at at 
least two different flow rate values of one fluid (blood or the 
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product of ultrafiltration) . 

In a further embodiment of the invention, the monitoring 
procedure includes determining the. values assumed by the blood 
characteristic in at least one zone of the blood circulation 
5 path and at at least two different flow rate values of two 
fluids (blood and the . product of ultrafiltration). 

In a further embodiment of the invention, the monitoring 
comprises a measuring stage of a blood characteristic, in a zone 
of the extracorporeal circuit arranged downstream of the blood 

10 withdrawal zone, or in a zone arranged upstream of the blood 
return zone, or in both above zones. The monitoring includes 
determining the blood characteristic in the vascular access, in 
the withdrawal zone, or in the return zone, or in both zones, by- 
means of one or more mathematical models describing the 

15 variation of the said blood characteristic between the zones of 
withdrawal and return in the vascular access and the measuring 
zones in the extracorporeal circuit. The mathematical models can 
be, in particular, models descriptive of the variation of the 
said blood characteristic in the passage through the arterial or 

20 venous needles. In an embodiment of the invention, these 
mathematical models comprise at least one parameter which is 
relative to the blood flow, or at least one parameter relative 
to the hematocrit of the blood, or both said parameters. In 
particular the mathematical models can be represented by one or 

25 more interpolating formulas of experimental data; the formulas 
can be, for example, second-order polynomials with one or more 
parameters chosen between the flow rate and the hematocrit of 
the blood. 

In a special operation of the invention, at regular time 
30 intervals the monitoring procedure determines the values assumed 
by the blood characteristic in at least one zone of the blood 
circulation path during the flow rate change, evaluates the 
variation of the blood characteristic, selects the values 
assumed by the blood characteristic when the variation has 
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exceeded a threshold limit value, and uses the selected values 
to calculate the value of: the blood characteristic at the 
vascular access. 

In a further special operation of the invention, at regular time 
5 intervals the monitoring procedure determines the values of the 
blood characteristic in two different zones of the blood 
circulation path during flow rate change, compares the variation 
of the blood characteristic detected in a first zone. . of the 
blood circulation path and the variation of the characteristic 
10 detected in a second zone thereof, selects the values of the 
blood characteristic when the difference between the variations 
has exceeded a threshold limit value, and uses the selected 
values in calculating the value of the blood characteristic at 
the vascular access. 

15 In another characteristic of the invention, in calculating the 
value of the characteristic of the vascular access, the 
monitoring procedure considers the values of the blood 
characteristic in a stationary blood flow situation, i.e. after 
having kept the flow rate constant for a determined period of 

20 time. 

The monitoring procedure is applied by means of a machine for 
blood treatment in an extracorporeal circuit, in particular for 
a machine for treatment of kidney failure, predisposed to 
perform one or more of the following therapies: hemodialysis, 
25 hemof iltration, hemodiaf iltration, pure ultrafiltration, 
plasmapheresis. 

The machine is provided with a timer for carrying out the 
monitoring procedure at least once during the extracorporeal 
treatment . 

30 The monitoring procedure can be initiated on command of an 
operator, or automatically at a predetermined moment during the 
treatment . 



- 16 / 74 - 



The extracorporeal circuit can be included in the complex of 
fluid distribution lines, of the disposable type, normally 
removably associated and used in a machine for treatment of 
renal failure. 

5 The machine is normally equipped with pressure transducers 
operating in the blood withdrawal line, before the blood pump, 
and in the blood return line, after the blood treatment unit. 

Further characteristics and advantages of the present invention 
will better emerge from the detailed description that follows, 
10 of a specific embodiment of the invention, illustrated purely in 
the form of a non- limiting example in the figures of the 
drawings . 

Brief Description of the Drawings . 

The description will be made herein below with reference to the 
15 appended figures of the drawings, here given by way of non- 
limiting illustration, in which: 

- figure 1 is a diagram of a machine for an extracorporeal 
blood treatment provided with a monitoring device of the 
vascular access according to the invention; 

20 - figure 2 is a diagram of blood flow in a patient connected 
up to the machine of figure 1 ; 

figure 3 is an electrical diagram which describes by 
analogy the circulation of extracorporeal and 
intracorporeal blood of the patient subjected to the 
25 extracorporeal treatment with the machine of figure 1 ; 

figure 4 shows a diagram of the relation between AP f and q hf 
where AP f = P af - P vf (difference between arterial pressure 
in the vascular access P af and venous pressure in the 
vascular access P vf ) and is the extracorporeal flow rate 
30 of the blood; 
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figure 5 is a diagram of the relation between (P vf - P v ) e 
q uf/ where (P vf - P v ) is the difference between the venous 
pressure in the vascular access P vf and the systemic venous 
pressure P v , and q uf is the ultrafiltration flow rate. 

5 Detailed Description. 

The machine illustrated in figure 1 is a machine for 
hemodiaf iltration comprising a unit for an extracorporeal blood 
treatment (a filter for hemodiaf iltration 1) having two chambers 
2, 3 separated by a semipermeable membrane 4. A first chamber 2 

10 has an inlet which is connected to an arterial line 5 (blood 
withdrawal line from the patient) of an extracorporeal blood 
circuit. The arterial line 5 is connectable with a vascular 
access 6 of a patient by means of an access tool constituted in 
the example by an arterial needle N A . The arterial line 5 is 

15 provided with a pressure sensor 8 and a positive displacement 
pump 9 for blood circulation along the extracorporeal circuit in 
the direction of the arrow 7. 

The first chamber 2 has an outlet connected to a venous line 10 
(blood return line to the patient) of the extracorporeal blood 
20 circuit. The venous line 10 is connectable to the vascular 
access 6 of the patient by means of an access tool constituted 
in the illustrated embodiment by a venous needle N v . The venous 
line 10 is provided with a pressure sensor 12. 

The second chamber 3 of the filter 1 ^has an inlet connected to a 
25 supply line 14 of a fresh treatment fluid (dialysis liquid) and 
an outlet connected to a discharge line 15 of a discharge fluid 
(the dialysis liquid and the ultraf iltered liquid) . The supply 
line 14 is provided with a supply pump 13 of the fresh treatment 
fluid. The discharge line 15 is provided with a drainage pump 16 
30 for the circulation of the discharge fluid in the direction of 
the arrow 11 . 

The dialysis machine further comprises a control and calculation 
unit 17 connected to a screen and also to a keyboard through 
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which the user communicates to the control and calculation unit 
the setting values for machine operation. One of the setting 
values which the control and calculation unit 17 receives, from 
the user is the blood flow rate in the arterial- blood 

5 withdrawal line 5. The control and calculation unit 17 can 
control the speed of the blood pump 9 in order to have the 
predetermined value of flow rate c^. The control and calculation 
unit 17 can be connected to at least one measuring device, able 
to provide information relating to the effective blood flow rate 
10 in the arterial line. The measuring device can comprise, for 
example, a flowmeter, or an encoder connected to the rotor of a 
blood pump. The control and calculation unit 17 is further 
connected to the pressure sensors 8 and 12 and receives 
therefrom the signals indicating the detected pressure. 

15 The control and calculation unit 17 controls the operation of 
the various motor devices of the machine, in particular the 
blood pump 9 and drainage pump 16, according to the instructions 
received from the user and the programmed algorithms contained 
in its memory. 

20 The machine can further comprise sensors (of known type and not 
illustrated) for detecting the blood viscosity upstream and 
downstream of the treatment unit 1. The sensors can comprise, 
for example, measuring devices for the blood hematocrit level. 

The control and calculation unit is programmed to carry out, 
25 automatically or by request of the user, a series of operations 
which enable the vascular access to be monitored. 

Figure 2 shows the patient's blood circulation subjected to 
extracorporeal treatment with the machine of figure 1. The 
vascular access 6, through which the extracorporeal blood 
30 circuit is connected to the cardio-vascular circuit of the 
patient is, in the embodiment, a fistula of the Cimino-Brescia 
type. In figure 2 H indicates the patient's heart, P denotes the 
pulmonary circuit, V denotes the vascular system (or systemic 
circuit, or intravascular circuit or intracorporeal circuit). 
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The arterial line 5 and the venous line 10 are connected at one 
end to the vascular access 6 and at the other end to the 
dialysis filter 1. 

Figure 3 shows an electrical diagram which, by analogy, 
describes the blood circulation of the patient subjected to the 
extracorporeal blood treatment. 

The legend to figure 3 is as follows. 

Quantities controlled by the control unit 17: 

Ck blood pump flow rate [ml/min] 

q uf ultrafiltration flow rate [ml/min] 

Known quantities (measurable directly or indirectly or 
determinable from indirect measurements using a mathematical 
model) : 

P am extracorporeal arterial pressure [mmHg] 
P™ extracorporeal venous pressure [mmHg] 

E art hydrostatic pressure related to the height level difference 
between the pressure sensor 8 in the arterial line of the 
extracorporeal circuit and the arterial needle N A [mmHg] 

E ven hydrostatic pressure related to the height level difference 
between the pressure sensor 12 in the venous line of the 
extracorporeal circuit and the venous needle N v [mmHg] 

R am hydraulic resistance of the extracorporeal arterial line 
[mmHg-min/ml] 

R™ hydraulic resistance of the extracorporeal venous line 
[mmHg • mi n /ml] 

P af vascular access arterial pressure [mmHg] 

P vf vascular access venous pressure [mmHg] 
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P a mean systemic arterial pressure (MAP) [mmHg] 

P v venous pressure (venous return pressure) [mmHg] . 

Unknown quantities to be determined: 

q a blood flow rate at the vascular access, upstream of the 
arterial access [ml/min] 

q f blood flow rate of artero- venous anastomosis in the 
vascular access tract comprised between the arterial access 
and the venous access, (q f = q a - q b ) [ml/min] 

q v blood flow rate downstream of the venous access, (q v = q a - 
q uf ) [ml/min] 

R d hydraulic resistance upstream of the vascular access 
[mmHg-min/ml] 

R f hydraulic resistance between the arterial access and the 
venous access [mmHg-min/ml] 

hydraulic resistance downstream of the vascular access 
[mmHg • min/ml ] 

In the diagram of figure 3 the extracorporeal blood circuit is 
traced in bold line, while the intracorporeal circulation in the 
vascular access is drawn in thin line. 

The nodes where the extracorporeal circuit meets with the 
vascular access are the zones where pressures P af e P vf are 
determined (either directly measured or calculated) . 

Various methods are known, based on mathematical models, for 
calculating pressures P af e P vf from known pressures P am e P^ in 
the extracorporeal circuit. Some of these methods are described 
in the scientific publications cited in the present description. 
Herein below details will be given of a method founded on a new 
mathematical model based on the "electrical diagram represented 
in figure 3 . 
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In the following a mathematical model is shown, also based on 
the electrical diagram of figure 3, representative of the 
hemodynamics of the vascular access of an extracorporeal blood 
circuit in which' the blood is removed from the patient through 
an arterial needle, is made to circulate through the 
extracorporeal circuit and is returned through a venous needle. 

The mathematical model describes the variation of pressure in 
the vascular access as a function of. the blood flow rate. 

The mathematical model is expressed in the following three 
equations which can be derived from the electrical diagram 
represented in figure 3 . 

a = _f a J- 

PafP,f-Rf\qa-q b ) 
p V f-P v = R v\9 a -g U f) 

where, as mentioned herein above, the symbols have the following 
meanings : 

q a = blood flow rate at the vascular access 6 (fistula) , 
upstream of the withdrawal point of the arterial needle N A 

qb = blood flow rate in the arterial line 5 of the 
extracorporeal circuit 

P a = mean systemic arterial pressure measured at patient's arm 

P af = arterial pressure in the vascular access 6, i.e. the 
pressure in the vascular access (in the embodiment, with a 
Cimino-Brescia fistula, this is a tract of arterialized 
vein) at the point of withdrawal of the arterial needle N A 

R d = resistance of the tract of arterialised vein comprised 
between the anastomosis and the point of withdrawal of the 
arterial needle N A 
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P vf = venous pressure in the vascular access 6, i.e. the 
pressure in the fistula at the return point of the venous 
needle N v 

R f = vascular resistance of the tract of fistula comprised 
between the two needles N A and N v and representing the 
resistance between the two points at which P af and P vf are 
determined 

P v = venous pressure of the blood in the distal venous branch; 
the Pv value can be unknown during the extracorporeal 
treatment; in this case it can be placed at a constant 
physiological value (e.g. P v = 0) 

1^ = vascular resistance in the venous branch of the blood 
return zone at the zone where venous pressure P v is 
evaluated; where P v = 0 , the resistance represents total 
venous resistance, i.e. the vascular resistance met by the 
blood in returning from the venous needle N v to the heart H, 
which constitutes an indicative value of the drainage 
efficiency of the venous circulation 

q uf = ultrafiltration flow rate (in case of hemodiaf iltration, 
q uf is the difference between the discharge fluid flow rate 
in the discharge line 15 and the fresh dialysis fluid flow 
rate in the supply line 14) . 

The pressures in the above -indicated mathematical model relate 
to atmospheric pressure. The arterial and venous pressures P af 
and P vf in the vascular access are measurable directly, for 
example using pressure sensors operating directly on the 
vascular access 6 in proximity or internally of the arterial and 
venous needles N A e N v . 

As previously mentioned, the pressures P af and P vf are also 
determinable indirectly using a mathematical model which 
includes, among its parameters, pressures P am and (arterial 
and venous pressures) measured in the extracorporeal circuit by 
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the pressure sensors 8 and 12 . The prior art comprises various 
mathematical models usable for calculating pressures P af and P vf 
when pressures P am and P^ are known. Some of the above -cited 
prior art contains examples of so-usable mathematical models. 
There follows a further example of a mathematical model usable 
for determining the intravascular pressures of the blood 
starting from the easily-measurable values of the extracorporeal 
blood pressures. 

Determination of and with and known. 

The mathematical model used comprises the two equations which 
can be derived from the electrical diagram of figure 3 : 

P af = P am + E art + R am ' <lb 

P vf = P vm + E ven ~ Km ' b ~ 9«f) 

Resistances R am and can be considered equal, with 

satisfactory approximation, to the hydraulic resistance of the 
arterial needle N A and, respectively, the venous needle N v ; it is 
therefore assumed for the sake of simplicity that the whole drop 
in pressure in the arterial and venous lines is concentrated at 
the respective needles. 

To calculate the hydraulic resistance R of a needle, the 
following mathematical model is used: it makes use of an 
equation which connects the hydraulic resistance of the needle 
with the blood flow rate and the blood hematocrit. 



R = (A 2 -q\ + 4 -q b +B 2 -Hci 1 +B X -Hct + B 0 )'R 



Poiseuille 



where 

qb = blood flow rate 
Hct = blood hematocrit 



R 



Poiseuille 4 

n r 
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L = length of needle 

r = radius of internal section of the needle 

R poiseuiiie is tne theoretical hydraulic resistance calculated using 
the Hagen-Poiseuille law for a liquid with viscosity equal to 
5 one . 

A 2 , A lt B 2 , B x and B 0 are coefficients characteristic of each 
needle, the value being obtained by means of experimental 
preliminary laboratory testing, by measuring the fall of 
pressure through the needle with different blood and hematocrit 

10 flow rates. In experimental tests the flow rate was varied 
within a range from 0 to 500 ml/minute, while the hematocrit was 
varied within a range from 30 to 45%. The coefficients differ 
for a same needle according to blood flow direction, that is 
whether the needle is used as an arterial needle or as a venous 

15 needle. These preliminary in vitro tests serve to experimentally 
characterise the needles which will then be used for the 
extracorporeal blood treatment. The tests include simulation of 
the extracorporeal treatment (for example dialysis) using a 
machine for performing the treatment (for example a dialysis 

20 machine) with an extracorporeal circuit lacking the device for 
effecting the treatment (for example lacking a dialyzer filter) , 
causing bovine blood to circulate, exiting from a container and 
returning thereto. The blood is kept at a constant temperature 
of 3 7 °C. The blood hematocrit is measured. The machine and the 

25 circuit used in the tests can be the same as those illustrated 
in figure 1. 

At intervals of about 1 minute the blood pump flow rate is 
changed, starting from a zero flow rate q^ 0 = 0 ml/minute and 
increasing it by 50 ml/minute up to a maximum flow rate of 500 
30 ml/minute (q bl =50ml/min, ^=10 0ml /min, q bi =i • 50ml /min , 

q blo =500ml/min) . In general, the flow rate assumes N different 
values qti with i = 0, 1, 2, N (N > 3). 

At each interval pressures P ami and are measured using the 
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pressure sensors placed along the extracorporeal circuit. From 
each pressure value measured, P ami and we subtract the 

hydrostatic pressure due to the different blood level in the 
container with respect to the point of measurement of the 
5 pressure on the machine. From pressures P ami and P^ we can 
deduce the pressure falls of the corresponding needles AP ai and 
AP vi , with i = 0, 1, 2, N (N > 3) . 

The same operations are repeated, each time controlledly 
changing the value of the hematocrit in the bovine blood. The 
10 blood flow rate values ^ are the same each time, i.e. = q^, 
with i = 0, 1, 2, N. 

The hematocrit can be varied by dilution with physiological 
solution (in this case the hematocrit diminishes each time) . For 
each series of operations the value of the hematocrit is 
15 measured. Purely by way of an example, the operations can be 
performed with the following hematocrit values: about 44%, about 
42%', about 40%, about 38%, about 36%, about 34%, about 32 %. In 
general the value of the hematocrit Hct assumes M different Hctj 
values with j = 1, 2, M (with M > 2) . 

20 Thus, for each needle we obtain a number N* M of values AP aij and 
AP vij with i = 0, 1, 2, N (with N > 3) e j = 1, 2, M (with 

M > 2) . 

A processor calculates the hydraulic resistances of the needle, 
normalised with respect to the Poiseuilie resistance, for one of 
25 the hematocrit values (for example Hct = Hct x ) according to the 
equation: 

K ail~ — 

Qbi K Poiseuilie 

K vi\ — 

Qbi K Poiseuilie 

in which 
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R ail = resistance of the arterial needle at flow rate = and 
with hematocrit Hct = Kct 1 

Rvi! = resistance of the venous needle at flow rate = q bi and 
with hematocrit Hct = Hct x 

5 AP ail = pressure drop on the arterial needle at flow rate = 
and with hematocrit Hct = Hct x . 

A p vii = pressure drop on the venous needle at flow rate 0^= 
and with hematocrit Hct = Hct x . 

Hence we obtain two series of values R ail and of resistances 

10 (one arterial and the other venous) corresponding to a 
determined hematocrit value (in the example Hct = Hct^ , with i 
= 0, 1, 2,..., N, with N = number of times we determine AP ail and 
AP vil at different flow rates o^.. 

Each of the two series of values (R a and RJ is interpolated by 
15 the processor using a second-order polynomial: 

R = A 2 q 2 b + A x q h + b x 

and we thus obtain, for each type of needle, a pair of 
coefficients A 2 and A x for each flow direction (i.e. we obtain a 
pair of coefficients which characterise the arterial needle and 
20 a pair of coefficients which characterise the venous needle) . 
Coefficient b x depends on the blood hematocrit value. 

Coefficients B 2/ B 1 and B 0 are obtained as follows. 

Let us for a moment consider only one blood flow direction 
through the needle: for example, the arterial needle. 

25 The processor also calculates the resistances R a of the arterial 
needle for the other hematocrit values Hct = Hctj (j = 2,..., M) , 
at different blood flow rates q b =q bi (i = 0, 1, 2,..., N) , thus 
obtaining various series of values: 
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Poiseuille 



These values of R a are interpolated, for each hematocrit value 
Hct, according to the blood flow rate q^, using a second order 
polynomial : 



R ai3 = d 2 ql +A l q b +b 3 for Hct = Hct 3 

R aiM= A 2<ll+ A l<lb+ b M f ° r HGt = H Ct M 

with i = 0, 1, 2, N (with N > 3), in order to obtain a series 
10 of values bj (j = 1, 2,..., M) . 

In substance, exemplifying the above-mentioned process step by 
step, for j = 1 the processor interpolates values R ail (for Hct = 
Hct x ) according to the equation 



15 and thus determines b x . 

Then it interpolates values R ai2 for j = 2 (for Hct = Hct 2 ) 
following the equation 



and determines b 2 , and so on up until j = M, thus obtaining M 
20 values of bj . 

At this point the processor makes a further interpolation, using 
the values of bj according to the equation 



5 



'ail 



A 2 <ll + Aq b +b 2 for Hct = Hct 2 



R ai\ = A 2<ll+A<l b +b x 



R ai2= A 2<ll+ A X<l b +b 2 



b = B 2 Hct 2 + B x Hct + B 0 



and thus determines coefficients B. 



'2 / 



B x and B, 
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The same series of interpolations is effected using the data 
relating to the venous needle. 



Hereafter we report some examples of values of the coefficients 
A 2 , A x , B 2 , Bi and B 0 experimentally obtained. 

5 With a needle having the following characteristics: gauge = 15 
(internal diameter = 1.6mm), length = 28mm, the following is 
obtained: 

A 2 (arterial) = -0.00004, A^arterial) = 0.0351, 

B 2 (arterial) = 0.0192, B x (arterial) = -0.9398, 

10 B 0 (arterial) = 21.059, R PoiS euiiie = 0.022 

A 2 (venous) = -0.000026, A ± (venous) = 0.0266, B 2 (venous) = 0.0403, 
Bi (venous) = -2.2937, B 0 (venous) = 41.969, R Po i se uiiie = 0.022 

With a needle having the following characteristics: gauge = 16 
(internal diameter = 1.4mm), length = 33mm, the following is 
15 obtained: 

A 2 (arterial) = -0.00004375, A x (arterial) = 0.0309, 

B 2 (arterial) = 0.0081, B^arterial) = -0.3226, 

B 0 (arterial) = 8.3882, R Poise uiiie = 0.0442 

A 2 (venous) = -0 . 00002875 , A 1 (venous) = 0.0193, 

20 B 2 (venous) = 0.0037, B ± (venous) = 0.0487, 

B 0 (venous) = 1.4565, R PoiS euiiie = 0.0442. 

The control and calculation unit 17 memory is preloaded with the 
coefficient values A 2/ A 1# B 2 , B x and B 0 of the most commonly used 
needles (the memory contains two series of coefficients for each 
25 needle, one for each blood flow direction, i.e. a series 
relating to a needle's use as an arterial needle and as a venous 
needle) . The control and calculation unit 17 recognises the 
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needle used in the extracorporeal treatment time by time and 
consequently in the calculation of P af and P vf uses the 
coefficients relating to the needle being used. Recognition of 
the needle can be automatic (for example by means of an 
identification system associated to the needle) or can be user- 
guided. 

Thus a mathematical model is defined, usable by the control and 
calculation unit 17 for determining the pressures in the 
vascular access by measuring the pressure in the extracorporeal 
circuit . 

Herein below some operative methods are defined by means of 
which a processor in the control and calculation unit 17 of the 
machine can monitor the vascular access during an extracorporeal 
treatment . 

First monitoring procedure. 

In this first operative mode is varied at q uf = constant ( = 
0) , while P am and P^ are measured. 

The- operative mode is now described step by step. 

a. Determine values P afl and P vfl of the arterial pressure and, 
respectively, the venous pressure in the vascular access 
(fistula) at a known blood pump flow rate q^. 

b. Save and store values q^, P afl and P vfl in a memory. 

c. Change the blood pump flow rate to a known value q^. At the 
same time the ultrafiltration flow rate q uf is kept 
constant . 

d. Keep the blood pump flow rate at q b2 for a determined period 
of time (for example about ten seconds) to let the system 
become stable . 

e. Determine values P af2 and P vf2 of the arterial pressure and, 
respectively, of the venous pressure in the vascular access 
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(fistula) at blood pump flow rate q^ . 

f. Save and store values q^, P af2 and P vf2 . 

g. Steps c-f can be repeated for a desired number of times so 
as to save and store a series of values q^, P afi , P vfi , with i 

5 = 1, 2, 3, N, where N is an integer number greater than 

1. 

h. Calculate R f and q a using the values stored in the memory 
and the mathematical model expressed by the equation 

Pa f -P V f=Rf{q a -q b ) 

10 i. Save and store the values calculated for R f and q a . 

j . Calculate R^, using at least a part of the stored values and 
the mathematical model expressed by the equation 

P vf -P v =R v \q a -q uf ) 
k. Save and store the calculated value for R v . 

15 1. Calculate R d using at least a part of the stored values and 
the mathematical model expressed in the equation 

■ *„ 

where P a (mean systemic arterial pressure or MAP) is 
measured at the patient's arm in known ways and the 
20 measured value of P a is transmitted to the control and 

calculation unit 17. 

m. Save and store the value calculated for R d . 

The calculation of R f and q a in point h can be done in the 
following way. 

25 The stored values of q^, P afi and P vfi , with i = 1, 2, N (with 
N > 2), are introduced into the equation 
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so as to obtain a system of N equations with 2 unknowns - q a and 

AP n =R r {q a -g bl ) 
AP f2 =R r {q a -q b2 ) 

where AP fi = P afi - P vfi with i = 1, 2, N (N > 2) 

The unknown quantities o^ and R f can be determined by calculating 
the optimal solution of the above- indicated equation system. 

If N = 2 the system has an analytical solution. 

If N > 2 the two unknowns q a and R f can be determined using an 
optimisation algorithm. 

For example the processor calculates the two values, one q a and 
the other R f , for which the corresponding values of AP f 
calculated by the above -indicated system of equations are the 
closest to the AP fi values previously determined at point e. 

The following calculation procedure can be used. Using the 
values stored in memory, q^, P afi and P vfi/ by means of a 
mathematical interpolation algorithm previously stored in memory 
the processor determines a linear equation which approximates 
the relation between AP f and q^. Then the value of at AP f = 0 
is calculated, using the above - indicated linear equation. The 
value of qb at AP f = 0 is assumed to be equal to the flow rate q^ 
of the vascular access. The value of q a thus determined is 
stored in memory. Further, the processor calculates the value 
assumed by AP f at q^ = 0 , once more using the same linear 
equation. The value of AP f at ^ = 0 is assumed to be equal to 
the product of R f -q a . At this point, using the previously- stored 
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value of q a the value of R f can be calculated with a simple 
quotient . 

Graph AP f - ^ of figure 2 illustrates this mode of procedure. 
The points in figure 2 represent the determined values AP fi of 
5 AP f according to the blood pump flow rate q^. The straight line 
interpolating the various points is the graphic representation 
of the linear mathematical relation which connects AP f with o^. 

The interpolation method can be any known linear interpolation 
method. The straight line of interpolation intersects the 
10 horizontal axis (qj at q a and the vertical axis (AP f ) at R f -q a . 

Another way of calculating q a and R f is based on the description 
of the relation between and AP f using a non-linear 

mathematical relation (for example a polynomial of a degree 
greater than one) , derived by the processor with an 

15 interpolation method using the values stored in the memory q^, 
p afi e p vf±- After having derived this non-linear relation, the 
value assumed by q^ at AP f = 0 is assumed to be equal to the flow 
rate q a of the vascular access. The value of q a thus determined 
is stored in memory. Further, the processor calculates the value 

20 assumed by AP f at q^ = 0, using the above-cited non-linear 
equation as well. The value of AP f at = 0 is assumed to be 
equal to the product of R f -q a . At this point, using the 
previously- stored value of q a it is possible to calculate, by a 
simple division, the value of R f . This value represents, in the 

25 embodiment, the value of hydraulic resistance R f at point = 0 
(i.e. at zero blood flow rate in the extracorporeal circuit). 

At point c, the blood pump flow rate is varied from to q^ 

so that, in consequence of the change of flow rate - qbi/ the 

pressure difference AP f = P af - P vf varies significantly in 

30 absolute value and sufficiently to be appreciated (for example 

at least 2 mmHg) , i.e. so that 

|AP fl - AP f2 | > 2 mmHg, 
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where 



APfi = P af i - P V fi and 

APf 2 = ^af2 - ^vf2 

The same occurs for each flow rate change from q^ to Pibd+D • Tne 
5 values of q^ are selected so that the difference between the 
minimum value and the maximum value of does not exceed a 

predefined value (for example about 600 ml/min) in order that 
and R f can be considered as constant in the calculation with good 
approximation. 

10 At point c. the ultrafiltration flow rate q uf is kept constant = 
0 . 

At point j. the resistance Rv is calculated assuming = 0. The 
R v stored in memory can be one of the estimated R v± or the mean 
value of the estimated R^. 

15 *w=^-^ 

At point 1. the resistance R d stored in the memory can be one of 
the R di calculated with equation (1) or the mean value of the 
calculated R di . 

P — P 

K di 

20 Second monitoring procedure. 

In the second operative mode q uf is changed to ^ = constant (not 
zero) , while P am and P^ are measured. 

The operative mode is now described step by step. 

a. Determine values P afl and P vfl of the arterial pressure and, 
25 respectively, of the venous pressure in the vascular access 

(fistula) at a known ultrafiltration flow rate q ufl at a 
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predetermined blood pump flow rate 
Save, and store values q ufl/ P afl and P vfl . 

Change the ultrafiltration flow rate to a known value q uf2 . 
At the same time the blood pump flow rate is kept 

constant and equal to the initial flow rate of point a. . 

Keep the ultrafiltration pump flow rate at value q uf2 for a 
determined period of time (for example about ten seconds) 
to let the system become stable.. 

Determine values P af2 and P vf2 of the arterial pressure and, 
respectively, the venous pressure in the vascular access 
(fistula) at ultrafiltration flow rate q uf2 of the blood 
pump . 

Save and store values q uf2/ P af2 and P vf 2 . 

Steps c-f can be repeated for a desired number of times so 
as to save and store a series of values q ufi , P afi/ P V f±/ with 
i = 1, 2, 3, N, where N is an integer number greater 

than 1 . 

Calculate q a and Ry using the values stored in the memory 
and the mathematical model expressed in the equation 

P vf -P v =K\q a -q uf ) 
Save and store the values calculated for and q a . 

Calculate R f using at least a part of the stored values and 
the mathematical model expressed in the equation 

P af -P vf =Rf{q a -q b ) 
Save and store the calculated value for R f . 

Calculate .R d using at least a part of the stored values and 
the mathematical model expressed in the equation 
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p - p ' 

a =— 2L 

m. Save and store the value calculated for R d . 

At point c, the ultrafiltration flow rate is changed from q ufl 
to Q U f2 so that, in consequence of the change in flow rate q uf2 - 
5 <5ufi# tne difference of pressure AP vf = P vf - P v significantly 
varies in absolute terms sufficiently to be appreciated . (for 
example at least 3 mmHg) , i.e. so that 

|AP vfl - AP vf2 | ^ 3 mmHg, 

where 

10 AP vfl = P vfl - P v and 

AP v f 2 = Pvf 2 ~ P V 

The same can be said for each flow rate change from •• q ufi to 

C[uf (i+1) • 

At point c. the blood flow rate in the extracorporeal circuit qj> 
15 is kept constant at a known value which is not zero. 

At point h. the calculation of and q a is performed in the 
following way. 

The stored values of q ufi/ P afi and P vfi , with i = 1, 2, N (with 
N ^ 2), are introduced in the equation 

20 P vf -P v =R v \q a -quf) 

so as to obtain a system of N equations with 2 unknown 
quantities qa and R„ . 

Pvr\-P v = R v \q a -q«n) 

p vfi-P v = R v\<ia -q„f2) 

25 
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P vJN~ P v= R v' a - QuJN ) 



The unknown quantities q, and Ry can be determined by calculating 
the optimal solution of the above -indicated equation system. 

If N = 2 the system has an analytical solution. 

5 If N > 2 the two unknowns q a and can be determined using an 
optimization algorithm. 

A calculation procedure which can be used is the following. 
Using the values stored in memory, q ufi and P vfi , the processor 
determines, by means of a mathematical interpolation algorithm 

10 previously stored in memory, a linear equation which 
approximates the relation between AP vf and q uf where AP vf = P vf - 
P v . Then the value assumed by q uf at P vf - P v = 0 is calculated, 
using the above -indicated linear equation. The value of q uf at 
AP vf = 0 is assumed to be equal to the flow rate q a of the 

15 vascular access. The value of q a thus determined is stored in 
memory. Further, the processor calculates the value assumed by 
AP vf at q uf = 0, once more using the same linear equation. The 
value of AP vf at q uf = 0 is assumed to be equal to the product of 
Rv-q a . At this point, using the previously- stored value of qa the 

20 value of can be calculated by a simple division. 

The plot of AP vf as a function of q uf in figure 5 illustrates 
this mode of procedure. The points in figure 5 represent the 
determined values AP vfi = P vfi - P v of AP vf as functions of the 
ultrafiltration pump flow rate q uf . The straight line 
25 interpolating the various points is the graphic representation 
of the linear mathematical relation which connects AP f with q uf . 
The interpolation method can be any known linear interpolation 
method. The straight interpolating line intersects the 
horizontal axis q uf at q a and the vertical axis of AP vf at R^-qa. 

30 At point j . (determination of R f ) the following procedure is 
observed. 

For each of the estimated values of P afi and P vfi , a corresponding 
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value of R fi is calculated using the above- indicated equation, 
from which it is obtained: 

P - P 

„ r afx r vfi 
K fi = 

The R f value stored at point k. can be one of the calculated 
values for R fi or the mean value of the R fi values. 

At point 1. (determination of R d ) the following procedure is 
observed. 

For each of the estimated values of P afi , a corresponding value 
of R di is calculated using the above -indicated equation: 

„ _ P a ~ P afi 



The R d value stored at point 1 . can be one of the calculated 
values R di or the mean value of the R di values . 

Third monitoring procedure. 

The equations which define the mathematical model of the 
vascular access used previously: 

9a R d 



Pa f -P V f=Rf{q a -q b ) 
P vf -P v = K\<ia -<iuf) 



can be reformulated so as to evidence the dependence of P a£ and 
P vf on P a , Ck, q uf and P v through the unknown parameters P^, R £ and 
Rv. The reformulated equations are as follow: 

P af = — - — - — P a — - — q b 4« f + — P v 

af R d +R f +R v R d +R f +R v R d +R f +R v Huf R d +R f + R v 
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p = Hi P + / ! q U + i r_ p 

vf R d +R f +R v .*. R d +R f +R v Hb R d +R f + R v Huf R d +R f +R v 

These equations can be rewritten as reported herein below. 

Paf = CaO ■ Pa + ' 9t + ^ 2 ' <7 u/ + (l ~ C aQ )• P v 

P yf = ^vo • P a + C vl ■ q b + c v2 ■ q uf + (l - c v0 )• P„ 
in which: 

R f + R v R d R f R d R v 

Ca0 R d +R f +R v ° al R d +R f + R v ° a2 R d +R r +R v 

R v R f R v _ R v -{R d +R f ) 



C V 1 = 



R d +R f +R v R d +R f +R v R d +R f +R v 

The third operating mode (as the following fourth and fifth 
operating modes) calculates at least a part of the coefficients 

10 c a0/ c al/ c a2 and c v0/ c vl/ c v2 and from these derives R d/ R f and R^. 
The calculation of the. coefficients is done starting from one or 
more known values for each of the following quantities: P a/ c^, 
q uf , P v/ P af and P vf . The quantities P a , q^^, q uf , P v are known 
through measurement. The quantities P af and P vf are known by 

15 direct measurement of the pressures in the vascular access, or 
by a process of calculation starting from the measurement of the 
pressures in the machine P am and P^. 

As the number of coefficients c a0 , c al , c a2 , c v0/ c vl/ c v2 is greater 
than the number of the resistances R d , R f and R^ , there exists a 
20 multiplicity of relations between the coefficients and the 
resistances. In general, knowledge of three coefficients enables 
a determination of the resistances. 

In the third operating mode both flow rates qt, and q uf are varied 
and the arterial pressure in the machine P am is measured, from 
25 which arterial pressure in the vascular access P af is calculated. 

In a specific embodiment in a first stage the pressure P am at 
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flow rates q^ = 0 and q uf = 0 is measured; in a second stage 
pressure P am at flow rates q^ ^ 0 and q uf = 0 is measured; in a 
third stage pressure P am at flow rates * 0 and q uf ^ 0 is 
measured. 

More in general, at q uf = constant (for example = 0) is varied 
and P am is measured at different values of q^. Thereafter q uf at 
q^ = constant (for example ^0) is varied and P am measured at 
different values of q uf . 

In this third operating mode a mathematical model of . the 
vascular access is used which is represented by one equation 
only: 

p af = c a0 • P a + c al • q b + c a2 - g tff ■ + (l - c a0 y P v 

from which coefficients c a0/ c al , c a2 can be derived, which are 
sufficient by themselves for the calculation of the three 
resistances R d , R f/ . 

In this third operating mode at least one measurement is taken 
of the patient's arterial pressure P a . Further, distal venous 
pressure P v is assumed to be zero; for this reason the equation 
used is simplified as follows: 

The third operating mode is now described step by step. 

a. Determine values P af0 of the arterial pressure in the 
vascular access (fistula) and the systemic arterial 
pressure of the patient P a0 at a known ultrafiltration flow 
rate q ufl = 0 at a predetermined blood pump flow rate qb = 0 . 

b. Save and store values P a0 and P af0 . 

c. Calculate c a0 by means of the equation 
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Save and store value c, 



Change the blood flow rate to a known value q^. At the 
same time the ultrafiltration flow rate q uf is kept constant 
and equal to the flow rate at point a. (= 0) . 

Determine values P afl and P al of the arterial pressure in the 
vascular access (fistula) and, respectively, of the patient 
at blood pump flow rate q^. 

Save and store values q^, P afl and P al . 

Steps d-f can be repeated for a desired number of times so 
as to save and store a series of values q^, P afi , P ai , with i 
= 1, 2, 3, N, where N is an integer number greater than 
or equal to 1 . 

Determine c al by solving the system of equations: 

P a fi ~ C a0 * Pai = C a\ ' Q bi 

with i = 1, N (N > 1) 

If N = 1 it is sufficient to solve a linear equation with 
only an unknown quantity. 

If N > 1 the value of c al is found by means of an 
optimisation algorithm which determines the optimal 
solution for the above -cited system. The searched- for value 
can be the value of c al which minimises the error between 
the values of P af calculated with the above system of 
equations, P afi *, where the asterisk * indicates that the 
value has been calculated, and the P afi values determined by 
measuring a pressure correlated with P af . The optimisation 
algorithm can be, for example, a linear regression 
algorithm. 

Save and store value c al . 

Change the ultrafiltration flow rate to a known value q ufl 
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not zero. At the same time the blood flow rate has a 
known value different to zero. 

Determine values P a£1 and P al of the arterial pressure in the 
vascular access (fistula) and, respectively, of the patient 
at ultrafiltration flow rate q ufl . 

Save and store values q^, q ufl , P afl and P al . 

Steps k-m can be repeated for a desired number of times in 
order to store, a series of values q uf j , P af j , P aj , with j = 1, 
2, M , where M is an integer number equal to or greater 

than 1 . 

Determine c a2 by solving the following system of equations 

Pqff ~ C a0 ' Paj ~~ C a\ ' 9 ' bk = C a2 ' 4 ujj 

with j = 1, 2,..., M (M > 1) 

If M = 1 it is sufficient to solve a linear equation with 
only an unknown quantity. 

If M > 1 the value of c a2 is found by means of an algorithm 
of optimisation which determines the optimal solution for 
the above system. The sought-after value can be the value 
of c a2 which minimises the error between the values of P af 
calculated using the system of equations P afj *, where the 
asterisk * indicates that the value is a calculated one, 
and the values of P afj determined through measuring a 
pressure correlated by P af . The optimisation algorithm can 
be, for example, a linear regression algorithm (as at point 
i . above ) . 

Save and store the determined value of c a2 . 

Determine R f , R^ and R d by solving the following system of 
equations which express the relat ion between c a0 , c al , c a2 and 
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R d =(l/c a0 -i)-(R f+ R v ) 



The value of the resistance R f can already be determined at 
step j . as both c a0 and c al are already known. 

r. Save and store the first determined values of R f , R v and R d . 

s. Determine q a using one of the equations of the mathematical 
model of the vascular access, for example: 

Pa" Paf 
a aj 

t. Save and store the value calculated for q a . 

In steps from 1. to n. the operation of measuring P aj can be 
omitted; in this case the values stored and used for the 
calculation are the same P ai values calculated at point h. at ^ 
= q^ and q uf = 0, or at point a. at q b = 0 and q uf = 0. 

Fourth monitoring procedure. 

Varying qb at q uf = constant (for example zero) and measuring P am 
and P^. 

In this case too we calculate at least a part of the 
coefficients c a0 , c al/ c a2 and c v0 , c vl , c v2 from which R d/ R f and 
are obtained. The calculation of the coefficients is done 
starting from the knowledge of one or more values for each of 
the following quantities: P a/ q^ t q uf/ P v/ P af and P vf . The 
quantities P a/ qb, q uf , P v are known by measurements. The 
quantities P af and P vf are known by direct measurement of the 
pressures in the vascular access, or by means of a calculation 
process which uses the measured values of pressures P am and P^ 

- 43 / 74 - 



* 



in the extracorporeal circuit . 

In the fourth operating mode the measures were taken at q uf = 0 
and we use a mathematical model which includes both equations of 
P af and P vf which in this case are simplified into the following . 
5 formulation: 

^V= c .p-^+^i-flf*+(i-^ 0 ).p F 

Pyf = Cvo * P a + c vl • q b + (l - c v0 ) • P v 

In the fourth operating mode the processor determines the four 
coefficients c a0/ c al , c v0 , and c vl and from these it calculates 
10 the three resistances R d/ R^, R f . 

In the fourth operating mode the pressures P af and P vf in the 
vascular access are determined, either by direct measuring or by 
measuring pressures P am e P^ in the extracorporeal circuit and 
calculating P af and P vf by means of a mathematical model. The 
15 pressures P af and P vf are determined at different values of the 
blood flow rate In the fourth operating mode, the arterial 

and venous pressures P a and P v of the patient are also considered 
in the calculation of the coefficients. 

As coefficients c a0 , c al , c v0 , and c vl are greater in number than 
20 resistances R d , R f and R v , there exists a multiplicity of 
relations between the coefficients and resistances. In general 
the knowledge of three coefficients enables determination of the 
resistances. It has been . found that the most precise 
determination of the resistances R d , R f and is obtained by 
25 using the three coefficients, c a0 , c al , and c v0 . 

The fourth operating mode is now described step by step. 

a. Determine pressures P af , P vf , P a , and P v with the blood pump 
flow rate and the ultrafiltration flow rate at nil (c^ .= 0 
and q uf = 0) . 

30 b. The values thus determined, P af0 , P vf0 , P a0 and P v0 , are stored 
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in memory. 

The processor calculates c a0 and c v0 by means of the 
equations: 




Change the blood flow rate, to a known value = ^ 0. 
Determine at least one value of P af , P vf# P a and P v when = 

Save and store values P afl/ P vfl# P al and P vl above -determined. 

Repeat steps from d. to f . for a predetermined number of 
times N in order to obtain a series of values q bi , P afi , P V fii 
P ai and P vi with i = l, 2,..., N (N > 1) . 

Calculate c al as a solution for .the system of equations 

P afi ~ C a0 ' P ai ~ ( l ~ C a0 )* P vi = C al ' 9ti 

If N = 1 the solution is immediate. If N > 1 the solution 
is obtainable with an optimization algorithm, such as for 
example a linear regression algorithm. 

Save and store the value of c al . 

Determine resistances R d , R f , and Ry by solving the 
following equations which express the relation between c a0/ 
c ai / c vo and R d , R f , Ry : 

R d =— c -^- 

C v0 C a0 

R f =-^— 
f c a0 -\ 
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^ C a\ ' C vO 

{CaO- C vo)'( C aO- 1 ) 

k. Save and store values R d/ R f , and Rv above - determined . 

1. Determine the flow rate of the vascular access q a using one 
of the equations of the mathematical model, for example the 
5 second: 

Paf-P*f=Rf(<la-<l b ) 

At point e., determination of the value of P v can be performed 
in two ways . 

The first consists in considering P v constant (P v = P v0 ) during 
10 variation in the blood flow rate q^, thus ignoring the 
variations in the venous pressure P v which actually occur during 
the various operative stages. Consequently the system of 
equations of point h. can be rewritten in the following way: 

Pqfi ~" C a0 ' Pai ~ 0- ~" C a0 )" ^vO ~ C a\ ' <J bi 

15 The second way consists in considering the variations in P v to 
be proportional to the variations in the arterial pressure P a/ 
thus : 




This is equivalent to assuming resistances R d/ R f and Rv to be 
20 constant during variation of q^ 

In this case the equation of point h. is: 

p 

P a fi ~ C aO 'Pai C ao) # ^vO t ~^~ C a\ ' Q ' bi 

Note that by substituting, in the above equation, c a0 with the 
expression 
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aO ~ 




vO 



vO 



as in point c. of the present operating mode, the following 
equation is obtained: 




aO 



C a\ ' <lbi 



5 which is the same equation that appears at point i. of the third 
operating mode, in which the contribution of P v was ignored. 

Fifth monitoring procedure. 

The fifth operating mode is similar to the third, with the 
difference that, instead of determining P af , P vf is determined. 

10 Briefly, the fifth operating mode consists in varying the blood 
flow rate while maintaining the ultrafiltration rate q uf 

constant, in varying the ultrafiltration rate while keeping the 
blood flow rate constant, and in determining the venous 

pressure in the vascular access P vf at various values of the 

15 above-mentioned flow rates. The processor determines the 
resistances R d , R f and 1^ and the flow rate q a in the vascular 
access by calculating the coefficients c v0/ c vl , c v2 using the 
equation 



20 and the operative stages cited for the third operative mode. 

The resistances are calculated by solving the following system 
of equations: 



P yf = C v0 ' P a + C v\ ' Qb + ^v2 * <luf + ( l " C v0 ) ' P v 



vO _ 



R d +R f +R v 



R d +R f +R v 
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R v \R d +R f ) 
Cv2 R d +R f +R V 

The flow rate of the vascular access q a is calculated as in the 
third operative mode. 

Note that, by means of the second monitoring procedure, and IV 
can be derived by determining two or more values for the venous 
pressure alone (P^ in the machine or P vf in the fistula) , with 
the equation 

P vf -P v = K'(<ia-quf) 

while for the calculation of the values of R f and R d< the values 
of arterial pressure (P am or P af ) are also used, as well as the 
other two equations of the mathematical model: 

Similarly a further monitoring procedure can be formulated on 
the basis of which the values of q a and R d are calculated, 
determining two or more values of only the arterial pressure (P^ 
in the machine or P af in the fistula), using the equation 

q " R d 

while for calculating the values of R f and the values of the 
venous pressure (P^ or P vf ) are also used, as well as the other 
two equations of the mathematical model : 

P vf ~P v =R v \q a -q uf ) and P af - P vf = R f \q a -q b ) . 

In all of the above -described modes, the measurements are taken 
with the system in a steady state. For example, the various 
measurements are taken after a certain time interval (for 
example about ten seconds) after the blood flow rate or the 
ultrafiltration rate has been changed. 
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Two numerical examples of the application of the invention are 
now reported. 

First example. 

This example uses the above -described first monitoring 
procedure, applied to the apparatus of figure 1. 

Direct measurement of pressures P a , P af/ P vf were taken at 
different flow rate values q^. The measurements taken are 



reported in the following table. 

<5b P a P af Pvf AP f 

(ml/min) (mmHg) (mmHg) (mmHg) (mmHg) 

300 100 51 42 9 

200 52 41 11 

100 54 40 14 

400 51 42 9 

500 50 43 7 



The equation of the straight line interpolating points AP f is as 
follows (see figure 4, where AP f is a function of : 

AP f = 0.016- (925 - q b ) 

From which the following values are calculated 
R f = 0.016 mmHg min/ml 
q a = 925 ml/min 

From the third equation of the mathematical model used (assuming 
P v = 0) we have qb! = 300 ml/min: 

p 

R =— — = 0.04 5 mmHg -min/ml 
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Given P a = 10 0 mmHg, for = 300 ml /rain we obtain: 
P - P 

R d =— ^1= 0.053 mmHg • mi n /ml 



Second example. 

The second example uses the fourth monitoring procedure. 

In the following the values of the pressure measured at 
different blood pump flow rates are reported. 

C[b ^af ^vf 

(ml/min) (mmHg) (mmHg) (mmHg) (mmHg) 
0 120 62 35 0 

150 118 59 37 

250 117 57 37 

350 114 53 38 

From these values we obtain: 

Pa/0 ~ PvO _ 

C~ = — = 0.52 

0 P _ P 

aO vO 
P —P 

By applying a linear regression algorithm to the following 
equation: 

P<xfi ~~ C a0 ' Pai ~ 0 C ao)' ^v0 = C a\ ' % bi 

the following value for coefficient c al was found: 
c al = -0.0155 

After which the following resistance values were found: 
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R d = 0.069 mmHg min/ml 
R f = 0.03 2 mmHg min/ml 
R v = 0.042 mmHg -min/ml 
From this we calculated: 

5 Q n ~ — — = 842 ml /mm 
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